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Abstract. We discuss the nonlinear rheology of dense colloidal dispersions at the glass transition. A first principles approach 
starting with interacting Brownian particles in given arbitrary homogeneous (incompressible) flow neglecting hydrodynamic 
interactions is sketched. It e.g. explains steady state flow curves for finite shear rates measured in dense suspensions of 
thermosensitive core-shell particles consisting of a polystyrene core and a crosslinked poly(N-isopropylacrylamide)(PNIPAM) 
shell. The exponents of simple and generalized Herschel Bulkley laws are computed for hard spheres. 
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INTRODUCTION 

Soft materials, such as particle dispersions, exhibit a wide range of rheological properties. While dilute colloids 
flow with a viscosity only slightly higher than that of the solvent, concentrated dispersions behave as soft amorphous 
solids. Under large deformations, they may yield and flow with a strongly shear-rate dependent viscosity, while their 
structure can remain amorphous and homogeneous |Q]] . 

In order to address the nonlinear rheology of such viscoelastic colloidal dispersions, recently, a first-principles 
constitutive equation under arbitrary time-dependent homogeneous flow was developed [2]. It includes transitions to 
arrested glassy states as function of the thermodynamic control parameters, thus allowing the interaction between slow 
structural relaxation and time-dependent external flow to be investigated. The approach generalizes the integrations 
through transients (ITT) approach developed for steadily sheared dispersions |3] to the general situation. It concen- 
trates on the thermodynamic stresses arising from potential and Brownian interactions between the particles, because 
they overwhelm hydrodynamic contributions close to vitrification and at low deformation rates, and because they are 
dominated by slow structural rearrangements, which can be approximated using mode coupling techniques |4|]. 

A colloidal model system of hard spheres was developed |5[], and experiments at high densities were performed 
10 to test the connection between glassy dynamics and nonlinear rheology, which is at the heart of the theoretical 
approach. To this end, linear frequency dependent shear moduli and nonlinear shear rate dependent flow curves were 
measured and analysed in parallel jg], Here, we summarize some of the observations on the steady state flow curves 
and provide some new insights into the flow curves close to the glass transition. 

THEORY 



The microscopic theoretical approach J2l starts with a system of /V spherical Brownian particles interacting via 
internal forces and dispersed in a solvent with a specified velocity profile v(r, t ) = K(t) ■ r. The time-dependent velocity 
gradient tensor K(t) is assumed spatially constant, thus excluding the non-homogeneous flows which occur in shear- 
banded and shear-localized states. While hydrodynamic interactions and fluctuations in the velocity profile are thus 
neglected from the outset, this Statistical Mechanics model can be investigated by Brownian dynamics simulations 02D, 
and constitutes a first microscopic approach to real glassy colloidal suspensions. The properties of this microscopic 
model have been worked out for low densities (partially even including hydrodynamic interactions [8]), and it provides 
the starting point for various (approximate) theories for intermediate concentrations. 

Considering, for the sake of discussion, the example of steady shearing [3], the aim of the microscopic approach is 
to compute the stationary many-body distribution function, which is possible albeit only approximately. It describes 
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FIGURE 1. (Color online) Flow curves of a suspension of colloidal hard spheres for two densities close to the glass transition, 
0eff = 0.580 (fluid state, diamonds) and e g = 0.608 (glassy state, circles) rescaled to collapse at high shear; from Ref. ((J. 
Parameters are T = 60D /R 2 H , y c = 0.1 and v a = 50k B T/R 3 H for <j) eS = 0.580, and T = S4D Q /R 2 H , y c = 0.105 and v a = llk B T IR\ 
for <^ e ff = 0.608, respectively, where Rjj is the hydrodynamic radius and Dq the single particle diffusion coefficient at infinite 

dilution; k B T is the thermal energy. Fits using the Fjj -model are given as dotted red lines [6]. The blue dashed cuves gives 
the generalized Herschel-Bulkley law Eq. Q} with m = 0.143, which holds asymptotically for yc y, at the glass transition; 
cr+ = OAOkfiT /Rff is the dynamic yield stress at the transition. The green solid curve is the Herschel-Bulkley law ((2J which holds 
for intermediate shear rates % <C y <C T at the transition; the exponent m = 0.473 follows for hard spheres using the Percus-Yevick 

approximation for the equilibrium structure factor |9, 4]. The magenta dot-dashed straight line labeled rjl denotes the (extrapolated) 
linear high-shear asymptote. The inset tests the exponent in, as da jdy x y 1_m <x const, should hold for the Herschel-Bulkley law. 



the complete probability distribution of the particle positions, and using it any stationary average can be computed. 
Its computation is achieved using formally exact projection operator manipulations and generalized Green-Kubo 
relations, followed by mode coupling approximations, which share similarities with mean field approximations y|]. 
The approximations implement the assumption that slow structural rearrangements dominate the thermodynamic 
stresses in dense dispersions, and that they can be approximated using (transient) density correlation functions. The 
later are taken from equations of motions closed with (again) mode coupling approximations. Competition of structural 
slowing down (called 'cage effect') with flow induced decorrelation ('flow advection' entering via Cauchy-Green 
tensors [2]) and cut-off of memory, determines the distorted microstructure and the stresses. Universal aspects of the 

(y) 

predicted rheology can be captured in 'schematic models', like the F^ -model, which retain the nontrivial dynamical 
scaling laws of the microscopic approach, but simplify its spatial structure. 



EXPERIMENTAL SYSTEM 

The particles consist of a solid core of poly(styrene) onto which a network of crosslinked poly(N- 
isopropylacrylamide) (PNIPAM) is affixed. Immersed in water the shell swells at low temperatures. Raising the 
temperature above 32°C leads to a collapse within the shell. Screening the remaining electrostatic interactions by 
adding salt, the system crystallises such that its phase diagram can be mapped onto the one of hard spheres[5]. The 
dependence of the packing fraction e ff (viz. the relative volume taken by the spheres) on the temperature is given by 
the hydrodynamic radius Rh determined from dynamic light scattering in the dilute regime. The flow curves were 
measured with a stress-controlled rotational rheometer MCR 301 (Anton Paar), equipped with a Searle system |g], 
upon decreasing the shear rate y and checking carefully for stationarity. 



COMPARISON OF THEORY AND EXPERIMENT 



The microscopic ITT equations contain a non-equilibrium transition between a fluid and a shear-molten glassy 
state. Close to the transition, (rather) universal predictions can be made about the non-linear dispersion rheology and 
the steady state properties. A qualitative change in the flow curves is predicted, where the stationary stress exhibits 
a Newtonian viscosity rj(y — ► 0) — > rjoy in fluid states, while it exhibits a dynamic yield stress a(y — > 0) — > <J + in 
glassy ones. Non-trivial power-laws in the flow curves exist close to the glass transition itself; for hard spheres, it is 
estimated to lie at e jf = 0.580 |6]. Right at the glass transition, a generalized Herschel-Bulkley law 

a(y^ 0) - a+ (1 + \7/%\ m + c 2 \Y/%\ 2 '" + \j/%? m ) (D 

holds and describes the flow curves over an appreciable part of the range y < %; the exponent is m = 0.151 for 
the exponent parameter A = 0.735 computed within mode coupling theory for hard spheres using the Percus-Yevick 
approximation. 

Beyond the asymptotic regime, a (simple) Herschel-Bulkley law appears in the flow curves of the F$ -model H 

o(t£%) = 6 Q + o l \i\ M . (2) 

Its exponent can again be calculated from X, which for hard spheres gives m — 0.473. This law provides a semi- 
quantitative fit of the flow curves for around two decades in / close to the glass transition as shown in Fig.[T] 

The dynamic yield stress <7 + can be read off by extrapolating the flow curve in the glass to vanishing shear rate. 
While this agreement supports the prediction of a dynamic yield stress in the ITT approach and demonstrates the 
usefulness of this concept, small deviations in the flow curve at low y are present. They indicate the existence of an 
additional slow dissipative process which may correspond to a melting of the glass at even longer times or smaller 
shear rates beyond the experimental window. 

CONCLUSIONS AND OUTLOOK 

Recently, the connection was emphasized between the physics of the glass transition and the rheology of dense 
colloidal dispersions, including in strong homogeneous flow |2|,|3|]. The equilibrium structure factor and one time 
scale (obtained from matching to short time motion and affected by hydrodynamic interactions) are required as inputs 
to the derived constitutive equation. Using model colloidal particles made of thermosensitive core-shell particles the 
predictions could be investigated in the vicinity of the transition between two stationary states under steady shear, a 
shear-thinning fluid and a shear- molten glass yllg]. 
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